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Abstract: The calculation of the magnetic moment of the deuteron in the framework 
\^ ', of the Bethe-Salpeter approach is performed. The relativistic corrections are calculated 

\ analytically and estimated numerically. It is shown that the main contributions are due to 

partial waves with positive energies and P-waves. A comparison with the non-relativistic 
schemes of calculations including mesonic exchange currents is made. 
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^ : 1 Introduction 

(N , 

The relativistic description of the deuteron has a long history. This problem still remains 
^ ! present although non-relativistic schemes of calculations are widely used to analyze reactions 

\ with the deuteron. There are different approaches which take into account relativistic effects. 

O I They can be separated into three groups: (i) the Bethe-Salpeter approach [|l|, (ii) a reduction 

Q ' of the Bethe-Salpeter equation (quasi-potential, light front dynamics), and (iii) prescriptional 

^ ■ treatments. 

The Bethe-Salpeter approach is the most consistent one. However, it contains some dif- 
ficulties that do not allow to establish a direct link to the non-relativistic calculations - e.g. 
the absence of a non-relativistic reduction of arbitrary kernels or the problem to interpret the 
^ ' abnormal parity states. 

C ■ Non-relativistic calculations with mesonic exchange currents give a reasonable description 

of the electromagnetic processes of the deuteron (e.g. elastic scattering and electro disintegra- 
^ \ tion). It is important to note that those mesonic exchange currents treated in a non-relativistic 

^ \ framework are partially included in the above mentioned relativistic calculations (i) and (ii) 

without additional assumptions. This has been demonstrated e.g. in ref. for the electro 
disintegration of the deuteron in the formalism of light front dynamics. 

Here, we present the investigation of a simple electromagnetic property of the deuteron, the 
magnetic moment, in the framework of the Bethe-Salpeter approach, and show the connection 
to the results of a non-relativistic calculation. 

2 Basic definitions 

The main object of the Bethe-Salpeter approach is the vertex function r(P, p) which obeys the 
Bethe-Salpeter equation, here written within the matrix formalism: 

[r(P,p)C]=./^t;(p,fc)r«^«(^ + A;)[r(P,A;)C]^(2)(^-A;)r(^), (1) 
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where C is charge conjugation matrix, v{p, k)V'^^^ ® F^^^ is the interaction kernel, and 5*^*^ S^^*) 
are Fermi propagators. The vertex function for the deuteron can be decomposed into basis 
states of definite values of orbital momentum L, spin S and p-spin p 0, 0]. The general form 
in the rest frame is given by 



(2) 



where are radial functions, and Tfj{p) are spin - angular momentum parts. 

In order to exhibit the p-spin dependence, the spin - angular momentum functions F^(p) 



may be replaced via FJj(p) = F^}^^''^(p), where 
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with a G {L, 5, J}, m the nucleon mass, p = 7^p'^, and F^^ given by 
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There are eight states in the deuteron channel (instead of two in the non-relativistic case), 
viz. 35++, 3^++^ 3^--^ 3^--^ 3pe^ 3po^ ipe^ 1 po (notation: 2-s+i^piP2). The normalization 

condition for this functions can be written as: 



P^ + P- 



P^ 
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+ Pw 
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+ Pipe + Pi po + Pspe + Pipe, 



(4) 



introducing pseudo-probabilities Pa that are negative for the states ^S^ , ^D^^ , ^Pf , ^P", ^Pf , 
^Pf, and positive for 'Si^, ^Di^ [|]. The calculation with realistic vertex functions give the 
following values: 



ipe _^3 po Ipe ^1 po 



4.8 -6-10-^ -0.88 -lO-^ -2.5-10-2 
5.1 -3.4-10-^ -9-10-2 -2.4-10-2 
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It is obvious that the main contribution to the normalization is due to the states with positive 
energies, and the contribution of the P-states is larger than that of the negative energies states 
by at least one order of magnitude. 
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3 The calculation of the magnetic moment 



The matrix element of the electromagnetic current used to evaluate the magnetic moment of 
the deuteron is given by 



{P'M'\J^\PM) 



le 



2tt^M 



(5) 



where \E'Af(P,p) = S^^\^ + p)Tm{P,p)S^'^\^ —p), the isoscalar form factors of the nucleons 
are given by Fl2{<f)i and p' = p + g/2, P' = P -\- q. The magnetic moment then is evaluated 
via 

1 -m _ /A/f' = -n \j \ M = n\ 

(6) 



Im /-,. (M' = +1|J,|M = 0) 
f^D = 7T V 2 lim , 

^ eM v^o ^y^TTrj 



with 7] = — g^/4M^. Since the integral vanishes at = 0, we expand the integrand of 
in powers of ^yfj. This is done in the Breit system. To this end the vertex functions need to be 
transformed from the rest system to the Breit system using the following formulae. 



(B) AB) ^ (B) 



where A(p(-^)) = (M + P^^'^-^o) / ^2M{E^b) ^ jvf), and the vectors P^^\ P' p^"', p are 
connected with the respective vectors in the rest system by 

The Lorentz transformation matrix C is of the form 
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From (0) it follows that p' = = + ig^) = £2^(0) + To simplify 

notation the vector will be denoted as = P = {po,Px,Py,Pz)- The components of the 
vector p' = p' are then given by 
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-- {l + 2r])p,-2^^1 + r]po + M^^l + r]. 
With the help of the transformations (|^), the integral (|^) may then be written as 



[P'M'\J^\PM) = ^ J d'pTrhM'{P^'\p')S^'^ 
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where the wave function \I/j\/(P(°\p) and the vertex function Tm'{P^'^\p') are taken in the 
deuteron rest frame. 

It is obvious from (|1^) that the sources of y/r] terms may appear in (i) the matrix A(P^^^), 
(ii) the propagator S^^\^P^^^ + p'), and (iii) the vertex function Tm'{P^^\p')- In detail this 
reads for the matrix (i) 



[A-'iP^""^^)]" = l + v^7o73, (li: 



and for the propagator (ii) 
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The corrections due to the vertex function (iii) may be separated into (a) corrections due 
to the radial wave functions, and (/3) due to the spin - angular momentum function. This will 
be demonstrated in the following using only one component of the vertex function (|^), i.e. 

TM'iP^'\p')=gip'oAp'\)^M'ip')- (13) 
For the radial part we get the following expansion 

9{Po, \p'\) = giPo, \P\) + VvP.i-^^^ + ^^^^^iUp^^ IpD- (14) 

The spin - angular momentum part can be written (e.g. for the ^Df^^ state) as 

rM'(p) = N{E'){m - p\)f{p\ ^)(m + p'^), (15) 



where 



P'l = Pi + v^(M-2po)[|7o-73], 

P2 = P2 + v^(M-2po)[|7o + 73], (16) 
- 2E^(iT^-2E(^ 



and corrections to r(p',^) can be calculated using (P^-(P^). 

4 Results of calculations 

The general formula for magnetic moment can be written as 

3 1 

/i+ = (/Xp + /i„)(P35++ + ^^3^++) - -{lip + f^n- 2)P^D++ + ^^ 
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Ail- = (A^P + l^n){P3p^ + Pipo) - -{PiPl + Pipo) - ^{P^Pl + P^P°) + ^1- 
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where Ra are the relativistic correction terms, viz 
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and iff'" (iff'" = iff) and G" are integrals of the form ^^j^^g / c^PoIpI with the integrands 
given by (compare with pi) 
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The relativistic corrections i?a can be estimated numerically for a separable model 0. They 
are smaller than the dominant terms by at least one order of magnitude. Taking into account 
different orders of pseudo-probabihties one gets 
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(1^ 



5 



3 1 

A/i = R+ + A/ii_ 
3 

Afii^ = --ilJ,p + fin) {Pip- + PiP° ) - (/ip + /in) (Pip^- + Pip° ) 

The contribution i?+ is negative and A/ii_ is positive. An estimation gives the following 
results 

R+/fiNR = -(5.8 - 11) X 10-2 %, A^i-/fiNR = 0.2 % (19) 

5 Conclusion 

We have shown that the expression for the magnetic moment in the Bethe-Salpeter approach 
can be written in a form closer to non-relativistic calculations. The additional terms in equation 
(|18D can be considered as relativistic corrections to the non-relativistic formula. 

The experimental value of the magnetic moment is known with a high accuracy: fiexp = 
0.857406(1). The non-relativistic value reflects only the D-state probability. Whereas in the 
relativistic corrections P-states play the dominant role. 

The magnitude of the corrections can be compared with the contributions of mesonic ex- 
change currents to the magnetic moment as extracted from ref. p|. The main contribution 
is due to the pair term, which leads to Afi/fi^R = 0.21 — 0.22% for different forms of the 



Bonn potential. The same size of this correction as compared to (|T^) may be considered as an 
indication that both corrections are of the same physical origin. 
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